The Interplay among Black Holes, Stars and ISM 
in Galactic Nuclei 

Proceedings IA U Symposium No. 222, 2004 

Th. Storchi Bergmann, L.C. Ho & H.R. Schmitt, eds. 



© 2004 International Astronomical Union 
DOI: 00.0000/X000000000000000X 



Long-Term Profile Variability of 
Double-Peaked Emission Lines in AGNs 

K. T. Lewis 1 , M. Eracleous 1 , J. P. Halpern 2 & T. Storchi-Bergmann 3 

1 Department of Astronomy and Astrophysics, The Pennsylvania State University, 525 Davey 
Laboratory, University Park, PA 16802, USA 

2 Department of Astronomy, Columbia University, 550 West 120th Street, New York, NY 

10027, USA 

3 Instituto de Fsica, UFRGS, 91501-970 Porto Alegre, RS, Brazi 

Abstract. An increasing number of AGNs exhibit broad, double-peaked Balmer emission lines, 
which arise from the outer regions of the accretion disk which fuels the AGN. The line profiles 
vary on timescales of 5-10 years. Our group has monitored a set of 20 double-peaked emitters 
for the past 8 years (longer for some objects). Here we describe a project to characterize the 
variability patterns of the double-peaked Ha line profiles and compare with those of two simple 
models: a circular disk with a spiral arm and an elliptical disk. 



1. Introduction 

Approximately 20% of Broad-Line Radio Galaxies (BLRGs) exhibit broad, double- 
peaked Balmer emission lines (Eracleous & Halpern 2003) . Double-peaked emission lines 
have also suddenly appeared in LINERs, most notably NGC 1097 (Storchi-Bergmann et 
al. 1993). The lines are thought to originate in the outer portions of the accretion disk 
(R ~ 100-10,000 T g , where r g = GM./c 2 ). The profiles of the double-peaked emission 
lines are observed to vary on timescales of 5-10 yr and are a manifestation of physical 
changes in the outer disk. Thus, the profile variations of these lines can be used to test 
various models for dynamical phenomena in the outer accretion disks of AGNs. 

Here, we compare the data of PKS 0921-213 with two simple, but physically motivated, 
models: a circular disk with a single precessing spiral arm and a precessing elliptical 
disk. Spiral arms are present in other astrophysical disks and provide a mechanism for 
removing angular momentum from the disk. An elliptical disk could result from a tidal 
perturbation of the disk by a massive object (e.g., a second black hole) or from the tidal 
disruption of a star. 

Spiral Arm: The spiral arm is modeled as a circular disk with a perturbation in the 
emissivity pattern. The disk has the following parameters: the inner and outer radii, £i 
and £ 2 (hi units of GAf./c 2 ); the inclination angle i; a broadening parameter er; and the 
slope of the axisymmetric emissivity law q (e oc The spiral arm perturbation is 

parameterized by: an amplitude A] pitch angle p; width 5; the inner radius of the arm 
£ sp ; and phase angle 4>. Elliptical Disk: The elliptical disk model is described by many 
similar parameters: i; <t; q; <j>; £i, and £2, where £ is the pericenter distance. The disk is 
circular up to £ e , then the eccentricity increases linearly to a value of e at £2- 

2. Characterizing the Profile Variability 

Both models involve many parameters and it is difficult to determine which region of 
parameter space will be most successful in reproducing the observed profile variations. 
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Figure 1. Variability of data and model profile properties with time. (Left) PKS 0921—213: 
The FWHM and FWQM, and corresponding shifts, are denoted with triangles and filled circles, 
respectively. (Center) Spiral Arm Model: Variation of profile properties with phase for a 
disk with a one-armed spiral. The model parameters are: q=1.5; i—50°; £=1200-8000r s ; a— 600 
km/s; A=5; p=15°; 5=35°; and £ sp = 1200r s . The FWHM and FWQM are denoted with solid 
and dashed lines, respectively. (Right) Elliptical Model: Same as the Center panel, but for an 
elliptical disk with q=1.5, i=50°, £=1200-8000 r g , a=600 km/s, e=0.15, and £ e = 1200r s . 



Thus, as a first step, we have begun characterizing the profile variability of our objects in 
a model independent way. Each profile is reduced to a set of easily measured quantities: 
the velocities of the red and blue peaks; the blue-to-red peak flux ratio; the full widths at 
half and quarter maximum (FWHM and FWQM); and the velocity shifts of the FWHM 
and FWQM centroids. Sets of model profiles, with a variety of input parameters, can be 
characterized in the same way for comparison with the data. 

As an example, in Fig. 1 we show the profile parameter variations for PKS 0921- 
213 (left), a spiral arm model (center) and an elliptical disk model (right). The circular 
disk parameters (£1,^2, 9,2, and a) for both models are tuned to best-fit parameters of 
the average profile of PKS 0931-213, while the additional model parameters are simply 
chosen for illustration. The most striking difference between the two models is that the 
elliptical disk profiles always vary smoothly with phase, while those of the spiral arm 
model vary much more sharply. However, both models show potential for reproducing 
the observed profile variation. By creating "libraries" of model profiles, we will be able 
to quickly select a small region of parameter space to test. 
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